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We present a search for Higgs bosons in multilepton final states in pp collisions at y/s — 1.96 TeV 
recorded with the DO detector at the Fermilab Tevatron Collider, using the full Run II data set with 
integrated luminosities of up to 9.7 fb~^. The multilepton states considered are eefi, e^f^, j-iThTh 
and like-charge e^/x^ pairs. These channels directly probe the HVV {V — W,Z) coupling of the 
Higgs boson in production and decay. The fj.ThTh channel is also sensitive to H ^ t^t~ decays. 
Upper limits at the 95% C.L on the rate of standard model Higgs boson production are derived in 
the mass range 100 < Mh < 200 GeV. The expected and observed limits are a factor of 6.3 and 8.4 
above the predicted standard model cross section at Mh ~ 125 GeV. We also interpret the data in 
a fermiophobic Higgs boson model. 

PACS numbers: 14.80. Ec,14.80.Bn 



I. INTRODUCTION 

The Higgs boson is predicted by the standard model 
(SM) as a consequence of the breaking of the electroweak 
symmetry, which gives mass to the weak gauge bosons. 
The ATLAS and CMS Collaborations at the CERN 
Large Hadron Collider have recently reported the obser- 
vation of a Higgs-like boson at a mass of Mh ~ 125 GeV, 
primarily in 77 and ZZ final states [l|, 0. Combining 
searches in the channel where the Higgs boson is pro- 
duced in association with a.W or Z boson, the CDF and 
DO Collaborations have found evidence for Higgs boson 
decay into bb pairs ^] . 

In this Article, we study final states with multiple 
charged leptons, including electrons, muons, and hadron- 
ically decaying tau leptons (t/j). We present the first 
Higgs boson search performed in the trilepton final states 
ee/i, e/i/i, and fiThTh with the DO detector. 

We also consider the production of like-charge e'^fi^ 
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pairs. This final state has the advantage of reduced back- 
ground from Z boson decay that is present in opposite- 
charge e+e", n'^iJ,~, and e^/z^ final states Q- This anal- 
ysis supersedes the previous searches in e^/i^ final states, 
which used integrated luminosities of up to 5.3 fb~^ 

The main Higgs boson production mechanisms relevant 
for this analysis are associated WH and ZH production 
and gluon-gluon fusion. The contribution from vector 
boson fusion is small and therefore neglected. The mul- 
tilepton channels are sensitive to Higgs boson decays into 
W^W^ and ZZ pairs, where the vector bosons {V) de- 
cay leptonically. These channels therefore directly probe 
the HVV coupling in production and decay. The trilep- 
ton searches are also sensitive to iJ — > r+r^ decays 
from associated production (WH, ZH) through hadronic 
tau decays in the iiThTh channel and through Icptonic 
tau decays in the ee/i and e/i/i channels. Searches for 
H T^T~ decays in final states with additional jets have 
also been performed using the full Run II data set Q. 

We also interpret the data in a fermiophobic model 
with a Higgs boson that does not couple to fermions 
and couples to W and Z bosons with SM strengths. 
Such searches have been conducted at the CERN e+e~ 
Collider LEP 0413 and at the Fermilab Tevatron Col- 
lider dH, The CMS Collaboration excludes fermio- 
phobic Higgs bosons with Mh < 124.5 GeV, 127 < 
Mh < 147.5 GeV, and 155 < Mh < 180 GeV at the 
95% C.L. in a model that assumes the couplings of the 
Higgs boson to other bosons are SM-like hj. 

Most results in this Article are based on the full Run II 
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data set collected with the DO detector at the Fermilab 
Tevatron and correspond to an integrated luminosity of 
9.7 fb^^. The analysis of the ^ThTh final state only uses 
data recorded after June 2006 with an integrated lumi- 
nosity of 8.6 fb~^. The results provide an important in- 
put for the combined Higgs boson search performed by 
the DO Collaboration [3] and for the Tevatron combina- 
tion [15]. 



II. DO DETECTOR 

The DO detector [l^ comprises tracking detectors and 
calorimeters. Silicon microstrip detectors and a scintillat- 
ing fiber tracker are used to reconstruct charged particle 
tracks within a 2 T solenoid. A liquid-argon and uranium 
calorimeter has a central section (CC) covering pseudo- 
rapidities \r]d\ up to « 1.1, and two end calorimeters 
(EC) that extend coverage to |77t;| w 4.2. The calorime- 
ters consist of electromagnetic (EM) and hadronic sec- 
tions segmented longitudinally in several layers. Muons 
are identified by combining tracks with patterns of hits in 
the muon spectrometer, which lies outside the calorime- 
ter and consists of a layer of tracking detectors and scin- 
tillation trigger counters in front of a 1.8 T toroid, fol- 
lowed by two similar layers after the toroid [31 . Trigger 
decisions are based on partial information from the track- 
ing detectors, calorimeters, and muon spectrometer. 



III. EVENT SIMULATION 

All background processes are simulated using Monte 
Carlo (MC) event generators, except the Zj back- 
ground in the e^/i channel and the multijet back- 
ground, which are determined from data. The Vl^-|-jets, 
Z/7* ^+^~-|-jet, and tt processes are generated us- 
ing ALPGEN with showering and hadronization pro- 
vided by PYTHIA [lO]. Diboson production {WW, WZ, 
and ZZ) and signal events are simulated using PYTHIA. 
All these simulations use the CTEQ6L1 [21] parton dis- 
tribution functions (PDFs). Associated production of 
Higgs bosons {WH and ZH) and gluon-gluon fusion are 
generated in 5 GeV increments of Mh in the range 
100 < Mh < 200 GeV. Tau lepton decays arc simulated 
with TAUOLA 22], which includes a full treatment of the 
tau polarization. 

Next-to-leading order (NLO) quantum chromodynam- 
ics calculations of cross sections are used to normalize 
the background contribution of ti [1^ and diboson 
processes. The WZ production cross section is corrected 
for Wj* interference using powheg 251- The cross sec- 
tion for W/Z+ jets production is normalized to a next-to- 
NLO (NNLO) calculation [2^. The transverse momen- 
tum (pt) spectrum of Z bosons is corrected to match the 
measured distributions The correction factor for the 
Pt spectrum of W bosons is the product of the Z boson 



correction factor and the ratio of the pt spectra of W 
and Z bosons calculated at NNLO [2§]. 

The cross sections for VH associated production are 
calculated at NNLO [2i,[33|. The NNLO calculation of 
Higgs boson production in gluon-gluon fusion takes into 
account resummation of soft gluons to next-to-next-to- 
leading- log (NNLL) [31]] . Higher order corrections to the 
Higgs boson production cross sections are computed with 
the MSTW 2008 PDF set [H. The simulated pr spec- 
trum of Higgs bosons from gluon-gluon fusion is corrected 
using the NNLO and NNLL calculation of hqt [H]. 
Branching fractions of the Higgs boson decays are cal- 
culated using HDECAY [11] . 

All MC samples are processed through a geant [bsI ] 
simulation of the detector. Data from random beam 
crossings are overlaid on the MC events to account for 
detector noise and additional pp interactions. The simu- 
lated distributions are corrected for differences between 
data and simulation in the reconstruction efficiencies and 
in the distribution of the longitudinal coordinate of the 
interaction point. 

To maximize signal acceptance, we use all events that 
pass our event selection without requiring a specific trig- 
ger condition. The residual efficiency loss from events 
that are not recorded depends on the event kinematics. 
We study ratios of kinematic distributions using the in- 
clusive trigger requirements and using a set of specific 
single lepton triggers. These ratios are then used to de- 
rive corrections on the normalization and shape of the 
kinematic distributions for MC events. 



IV. OBJECT IDENTIFICATION 

The signal comprises electrons, muons, and tau leptons 
that are isolated from other particles in the detector. 

Electrons are characterized by their interaction and the 
resulting shower shape in the EM calorimeter. The elec- 
tron clusters in the EM calorimeter are required to match 
a track in the central tracker. The energy is measured 
in the EM and the first hadronic layers of the cal orime- 
ter within a cone of radius 7^ = ^/{Arj)^ + (A0)2 = 0.2, 
where (j) is the azimuthal angle. The electron cluster must 
satisfy a set of criteria: (i) calorimeter isolation fraction, 
/iso = (Etot - Eem)/Eem, less than 0.15 for the CC 
region and less than 0.1 for EC, where Etot is the to- 
tal energy in the cone of radius TZ = 0.4 and Eem is 
the EM energy in a cone of radius TZ = 0.2; (ii) frac- 
tion of the EM energy to the total energy greater than 
0.9, and (iii) ratio of the electron's transverse momen- 
tum measured by the calorimeter and by the tracking 
detector, respectively, less than 8 (CC only). In addi- 
tion, the value of an eight- variable likelihood for electron 
candidates is required to be >Ce > 0.05 [1]. The electrons 
must also satisfy a requirement on a neural network dis- 
criminant with seven input variables in the CC and three 
in the EC region, including isolation and shower shape 
variables to improve the discrimination between jets and 
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electrons. The sum of the charged particle tracks' pt in 
an annulus of 0.05 < TZ < 0.4 around the electron direc- 
tion must be less than 3.5 GeV in the CC and less than 
(-2.5 X \r]d\ + 7) GeV in the EC. 

Muons are identified by the presence of at least one 
track segment reconstructed in the muon spectrometer 
which is spatially consistent with a track in the central 
detector, where the momentum and charge are measured 
by the curvature of this track. Muon isolation is im- 
posed with two isolation variables defined as the scalar 
sums of the transverse energy in the calorimeter in an 
annulus of radius 0.1 < TZ < 0.4 around the muon direc- 
tion and of the momenta of charged particle tracks within 
TZ — 0.5. Both variables, divided by the muon px, must 
be less than 0.2. To reduce the effects of charge mis- 
reconstruction, additional selection criteria on the track 
quality are applied in the e^p^ channel. 

Three types of tau lepton decays into hadrons are iden- 
tified by their signatures. Type-1 tau candidates consist 
of a single track and its associated energy deposit in the 
calorimeter, without any additional separate energy de- 
posits in the EM section. This signature corresponds 
mainly to — >■ n^iy decays and also includes leptonic 
e^yu decays. For type-2 tau candidates, we re- 
quire a track and its associated calorimeter energy de- 
posit, plus a separate energy deposit in the EM calorime- 
ter consistent with a tt*^ — 77 decay, as expected for 
TT^TT^v decays. Finally, type-3 tau candidates con- 
sist of two or three tracks, combined with an energy de- 
posit in the calorimeter. This corresponds mainly to the 
decays 7r^7r^7r^(7r°)i'. In this analysis, type-3 tau 

candidates are required to have three tracks and an asso- 
ciated net charge of ±1. For each tau-type, a neural net- 
work is designed to discriminate r^, from jets. The neural 
network discriminants are required to be NNr > 0.75 for 
types-1 and 2, and NNr > 0.95 for type-3 [36]. The input 
variables for these neural networks are based on isolation 
variables for objects and on the spatial distributions of 
showers. 

Variables that include information on the imbalance 
in transverse energy (T^t) caused by neutrinos are used 
to improve the discrimination between signal and back- 
ground. The is calculated using the transverse en- 
ergy measured in the calorimeter, corrected for the pres- 
ence of identified muons. Two modified l^x variables, 
and S{]^j^), are used to reject events where the T^j^ 
arises from detector effects and not from neutrinos. In 
events where the opening angle </> between the l^x direc- 
tion and the nearest lepton or jet is small, the resolu- 
tion of the l^rp measurement is dominated by the uncer- 
tainty on the measured lepton or jet energy. Less signif- 
icance is assigned to this region by using ^j', defined as 
$Jrp = sm(p a (j) < tt/2 and — elsewhere. The 
significance S{}/!rp) \37} is defined so that larger values 
of S{$j') correspond to measurements that are less 
likely to be caused by fluctuations in jet energies. 

Jet variables are used to discriminate between signal 
and background in the efip and e^/i^ channels but not 



in the event selection. We identify jets using a midpoint 
cone algorithm with a cone size of TZ = 0.5, based 
on energy deposits in the calorimeter. We require p^^* > 
15 GeV in the e/i/i channel and p^^* > 20 GeV in the 
e^//* channel. In both cases the jets must lie within 
1^?^ I < 2.4. 

V. EVENT SELECTION 

The event selection is designed to maximize sensitiv- 
ity to a SM Higgs boson signal in each channel sepa- 
rately. The leading muon {pi) in the e/x/x and ^ThTh 
channels, the leading electron (ei) in the ee/i channel, 
and the electron in the e^/i^ channel are required to 
have > 15 GeV. All other selected leptons must have 
a transverse momentum of px > 10 GeV. The pseudora- 
pidity of at least one of the selected muons in all channels 
except ee/i and of both Th candidates in the HTuTh chan- 
nel must be \'q^\ < 1.6 and \r]'^\ < 1.5, respectively. The 
transverse momentum of type-1 and type-2 candidates 
must be pr > 12.5 GeV, and we require px > 15 GeV 
for type-3 Th candidates. 

The leptons in the events originate from a pp inter- 
action vertex, which is required to have a longitudinal 
coordinate located within 60 cm of the nominal center of 
the detector. The maximum difference between the lon- 
gitudinal coordinate at the distance of closest approach 
to the beam axis for all lepton pairings in an event must 
be less than 3 cm. 

To facilitate combining channels, we ensure that there 
is no overlap between them. All events with at least two 
electrons and at least two muons (ee/z/x) are included 
in the ee/i sample and removed from the e/^/i sample. 
All events included in the other trilepton final states are 
removed from the pThTh sample. We also reject events 
with an additional electron or muon in the e*/x* channel. 

We construct a variable M{Ul]prp) that is the invariant 
mass of the three leptons and the , where the vec- 
tor is assumed to have a longitudinal momentum compo- 
nent equal to zero. We require M{M$rp) > 100 GeV for 
the ee/x and e/x/i final states to reject .Z-l-jets background. 
To remove Z — > pjj, events with final state radiation, we 
require e/i/i events in the range 75 < M^ejifj) < 105 GeV 
to have $t > 20 GeV. In the fiThTh final state, we require 
^T > 15 GeV and a transverse mass Mt{p) > 20 GeV. 
The transverse mass 

Mt{£) = y^2p^^T(l - cos?!))) 

is calculated using the azimuthal angle 4> between the 
charged lepton {£ = e, /i) and the direction of the 

Muons and electrons must be separated from jets by 
7^ > 0.1 in the epfi, ee/i, and e^/x^ channels. All selected 
leptons are required to be separated by 7?. > 0.3 from each 
other. This is increased to 7^ > 0.5 for the pairings of 
T/j candidates and for the pairing between candidates 
and muons. The sum of the charges in the pr^Th final 
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TABLE I: Numbers of events in data, predicted background, and expected signal for Mh = 125 GeV after the event selection. 
The numbers are shown for the different samples separately, together with their total (statistical and systematic) uncertainties. 
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FIG. 1: (color online). Distributions of (a) the product of the muon charges, g^j^ x g^jj , (b) the scalar sum of transverse momenta 
and -^y, "^Pt + (c) the transverse momentum of the electron, py, (d) the significance 5(.^y), (e) the maximal Mt2 for 
all lepton pairings, max{MT2 (-^^-^^'-^^t)}, and (f) the pseudorapidity of the trilepton system, 77'^^'^. The distributions are shown 
after the event selection. The e/ipA,B,c samples are shown in the left, middle, and right column, respectively. The data are 
compared to the sum of the expected background and to simulations of a Higgs boson signal for Mr ~ 125 GeV, multiplied 
by factor of 10 for the efifiA channel and 75 for the other channels. 



state must be ±1. The electron and muon are required to 
have the same charge in the e^/z^ final state. No lepton 
charge requirements are applied for the eefj, or efifj, final 
states to maximise the sensitivity to signal. 

We divide the e/x/x channel into three samples with dif- 
ferent signal and background composition to increase sen- 
sitivity to a Higgs boson signal. The e/i/iA sample con- 



tains events where the dimuon mass is outside the range 
60 < M{fi^) < 130 GeV and all events with like-charge 
muons. The second sample (e/i/XB) contains events with 
oppositely charged muons, 60 < M(/i/i) < 130 GeV, and 
S{lprp) > 2. The e/i/xc sample consists of all remaining 
events with 60 < M(/i/i) < 130 GeV and S{^t) < 2. 
The number of events in data, the expected back- 
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FIG. 2: (color online). Distributions of (a) trilepton transverse momentum, p^'^, and (d) the minimal transverse mass, 
min{MT(ei), MT(e2)}, for the eefi samples, (b) the transverse momentum of the leading th candidate, , and (e) the di-tau 
invariant mass, M{ti,T2), for the ^r^Th samples, and (c) missing transverse energy, and (f) the invariant mass of the 
electron and muon, Mep, for the e^/i* samples. The distributions are shown after the event selection. In addition, the final 
selection step requiring the output of the first BDT to be > 0.3 and min{i\fT(e), MT(At)} > 7 GeV has also been applied to the 
e^H^ sample. The data are compared to the sum of the expected background and to simulations of a Higgs boson signal for 
Mh = 125 GeV, multiphed by a factor of 10. 



ground and signal, after all selection criteria described 
in this Section have been applied, are given in Table ID 



VI. INSTRUMENTAL BACKGROUNDS 

Instrumental backgrounds are caused by leptons pro- 
duced inside jets, low-multiplicity jets that are recon- 
structed as Tfi candidates, photons or jets misidentified 
as electrons, and by opposite-charge efi pairs where one 
of the charges is incorrectly measured. 

The M^-l-jets background in the e/i/i and ee/i samples 
is expected to be small, and its contributions is there- 
fore described only by the simulation using the theoret- 
ical cross section. Since the W^-t-jets background is ex- 
pected to contribute more for the fiThTh and e^fi^ final 
states, their normalisation is obtained using data and 
then applied to the simulated kinematic distributions of 
the W-l-jets events. 

To model the VF-|-jets background for fiThTh final 
states, we select a data sample enriched in W^-|-jets 
events. We require that events pass all selection criteria, 
except the requirements on the NNt outputs. In addi- 
tion, we require MTi^J-) > 40 GeV, and plji > 20 GeV. 
The normalization factors that are applied to the simu- 
lation are determined from the ratio of the event yields 
in the W^-|-jets enriched region for data and simulation. 
They are determined separately for each type of Th candi- 
date, and for same-charge and opposite-charge pairs. 



To normalize the simulated M^-|-jets background in 
the e^/z^ final state, we select a data sample enriched 
in VF-|-jets events by requiring > 20 GeV, an 

inverted electron likelihood criterion £e < 0.7, and 
min{MT(e), Mt{h)} > 20 GeV. This data sample is used 
to derive separate normalization factors for jet multiplici- 
ties A'jct = 0, 1, and > 2 that are applied to the simulated 
M^-l-jets background samples. 

Multijet background is negligible for the e^/z, ee/i, and 
fiTiiTii final states. The multijet background in the e^jjL^ 
channel is determined using a sample with relaxed lep- 
ton identification requirements. The events in this sam- 
ple contain one electron with lower requirements on the 
likelihoods and neural networks compared to the stan- 
dard electron identification used in the event selection, 
and one muon without isolation requirements applied. 

To obtain the correct normalization of this sample, we 
calculate separate fake rates for electrons and muons, 
given by the ratio of the number of events with a stan- 
dard lepton to the number of events with a fake lepton. 
Fake leptons are defined by applying the relaxed lepton 
identification requirements but removing standard lep- 
tons. This fake rate is calculated separately for electrons 
and muons as a function of their transverse momentum 
in a sample with < 15 GeV. This selection ensures 
that both samples used to calculate the fake rates are en- 
riched in multijet events. The shape and normalization 
of the multijet contribution is then obtained by applying 
the product of the fake rates for electrons and muons to 
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the multijet data sample selected with the relaxed lepton 
identification requirements. 

Photons reconstructed as electrons contribute to the 
background in the e^/x final state. To estimate the con- 
tribution from the background, we select events with 
two muons, and with a photon or electron (labeled e/7) 
identified using identical criteria, apart from the track 
matching requirement, which is reversed for the photon. 
In addition, we require iflrp < 20 GeV. The invariant mass 
of the two muons and the e/7 has to be in the range be- 
tween 75 and 105 GeV. This sample is enriched in events 
with final state radiation from Z /^* — )• /i^/i^ decays. A 
normalization factor is calculated as the ratio of the num- 
ber of events with electrons to photons in this sample. To 
estimate the Z7 background, this normalization factor is 
applied to a 7/x^ data sample selected in the same way as 
the e/i/i sample, except for the reversed track matching 
requirement used to define the photon. 



VII. MULTIVARIATE DISCRIMINANTS 

Boosted decision trees (BDTs), as implemented in the 
TMVA package [HI , are used to discriminate between sig- 
nal and background. The BDTs are trained for each mass 
point separately in the range 100 < Mh < 200 GeV, in 
steps of 5 GeV. For each background or signal process, 
the event samples are split into subsamples for training 
the BDTs. The BDTs are then applied to the subsamples 
not used in the training to derive limits on the Higgs bo- 
son production cross section. The only exception is the 
background in the e/x/x final state, where both simu- 
lation and data are used. We train the BDT with a Z^ 
MC sample and then apply the BDT to the kinematic 
distributions estimated using the Z^ data sample to ob- 
tain the BDT distribution used in the limit setting. This 
procedure reduces fluctuations in the training caused by 
the small number of data events. 

The BDTs exploit kinematic differences between Higgs 
boson production for a given Mh and background. The 
variables used as inputs to the BDTs are given in Ta- 
ble [ni They are based on the transverse momenta of the 
leptons and jets, the ^j^, angular variables, charge corre- 
lations between leptons, and on the invariant masses of 
the pairings of leptons, jets, and the $rp- Jet variables 
are used to calculate the discriminants in the e/i^ and 
e^H^ channels only. The variable Mt2 is an extension of 
the transverse mass Mt to final states with two visible 
and two invisible particles [13, Slj . Some of the variables 
use constraints given by the Z boson mass. On average, 
the opening angle between leptons in H WW decays 
are smaller than for background and their direction is 
opposite to the direction of the l^x because of spin cor- 
relations in the decay of a scalar Higgs boson. Other 
variables, such as the likelihood Ce, reject events with 
misidentified leptons. Distributions of some of the in- 
put variables used for BDT training are shown in Figs. [T] 
andO 



A single BDT for each mass point is used to dis- 
criminate between signal and all background processes 
for the ee/i sample and for each of the three e/i/i sub- 
samples. The output distributions of the BDTs, shown 
in Fig. [3Ua)-(d) for data, signal with AIh — 125 GeV, 
and for expected background, are used to discriminate 
between signal and background. 

Two BDTs are trained for the fJ.ThTh sample, where 
the first BDT discriminates between signal and all back- 
ground sources except diboson production and the sec- 
ond BDT between signal and the dominant diboson back- 
ground. Events that pass a selection requirement on the 
first BDT discriminant of > 0.680-0.788, determined sep- 
arately for each Mh value to optimize the discrimination 
between signal and background, are used as input to the 
second BDT. The output distributions of the first BDT 
with Mh = 125 GeV is shown in Fig. [Hfe) and for the 
second BDT in Fig. [3Jf), using all events where the first 
BDT output is > 0.744. The output of the first BDT is 
used as the discriminant in the limit setting for all events 
that fail the requirement on the first BDT output and the 
output of the second BDT for all remaining events. 

The output distribution for the first BDT used for the 
e^/j,^ channel is shown in Fig. [3fg). It discriminates 
mainly between signal and M^-l-jets as well as multijet 
production. After requiring the output of the first BDT 
to be > 0.3 and mm{MT{e), MxilJ.)} > 7 GeV in a fi- 
nal selection step, the number of expected background 
events at Mh = 125 GeV is reduced from 809 ± 93 to 
122 zL 7, while the expected number of signal events is 
only reduced by 23% (see Table [HJ). A second BDT is 
trained to discriminate between signal and the remain- 
ing background sources, which are mainly from diboson, 
W-l-jets, and .Z-l-jets production, in the remaining events. 
The output of the second BDT for this sample, shown in 
Fig. [3jh) for Mh — 125 GeV, is used as discriminant in 
the Higgs boson searches. 



VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties on both background and sig- 
nal, including their correlations, are taken into account 
as uncertainties on the normalizations and on the shape 
of differential distributions. The theoretical uncertainty 
on the cross sections for Z/j* £'^£~, W+^ets, tt, and 
diboson production vary between 6% and 7%. The theo- 
retical uncertainty on the VF-|-jets cross section is applied 
for the ee/i and e/i/i channels, where this background is 
normalized using the theoretical prediction. The uncer- 
tainty from the normalization of the Vl^-fjets background 
in the in the fiThTh and e^/i^ channels is 5% — 6%. Ad- 
ditional shape dependent uncertainties on the H^-l-jets 
distributions are applied in the e^/x* channel to account 
for uncertainties in the description of the pt spectrum of 
W bosons as well as initial and final state radiation. The 
systematic uncertainty on the normalization of the mul- 
tijet background in the e^/i^ channel is estimated to be 
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TABLE II: Set of variables used in training of the BDT for each final state. The charges of the leptons are qi {£ = e,fj,). The 
angle (j){l\,l2) is taken between the two leptons, and the angle (j>{ll,l') between the dilepton system [U) and the lepton with 
the different flavour {I'). The variables TZ{£,£) and Mt2{U' IPx) calculated for all lepton pairings. The pairing with the 
smallest and largest values are denoted by min{ } and max{}, respectively, and mid{} corresponds to the third pairing. The 
mass, transverse momentum, and pseudorapidity of the trilepton system are denoted by M(£££), p^^, and rj^^^, respectively, 
and ^ is the scalar sum of the transverse momenta of the three leptons. The variable fcp is the fraction of charged particle 
tracks associated with the jet that point back to the same vertex as the leptons. If no jets are present, the jet variables are 
set to zero. All variables given for the jiThTh channel are used in the first BDT, whereas only pp, M{tt), Mri^y), and are 
used in the second BDT. All variables given for the e^/i* channel are used in both BDTs, except py* and ^(jet, which 
are used in the second BDT only. 
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20% by studying its dependence on jet multiplicity and 
lepton Pt- 

The theoretical uncertainty on signal production cross 
sections is 5% for gluon-gluon fusion and 6.2% for asso- 
ciated VH production. The uncertainty on the measured 
integrated luminosity is 6.1% ^^]. The systematic uncer- 
tainty on the lepton identification efficiency is 1.8% — 4% 
per muon and 2.5% per electron. The total uncertainty 
on the identification efficiency for both Th candidates, 
including the uncertainty on the neural network discrim- 
inant used to distinguish Th candidates from jets, is 7% 
per event. An uncertainty on the normalization of the 
signal and the simulated background is derived by com- 
paring distributions of data obtained using the inclusive 
trigger method with samples obtained using only single- 
lepton triggers. The resulting uncertainty is 1.5% — 5%. 

The uncertainty on the signal acceptance from the un- 
certainty on the parton distribution functions is 2.5%. 
The uncertainty on the probability that leptons originat- 
ing from jets are selected in the eefi and efifi channels is 



30% for the Z/W+iets and the WW samples. The un- 
certainty on charge misidentification is 20% for the e/i/i 
channel and 16% for the Z-|-jets background and 50% for 
the tt, and WW background in the e^^fj,^ channel. A 
7.3% systematic uncertainty is assessed on the Z7 back- 
ground in the e/z/z final state. The uncertainties on the 
description of the px distributions of W and Z boson, 
the uncertainties on the px resolutions for electrons and 
muons, and the uncertainties from jet energy resolution 
and efficiencies are found to have a negligible effect on 
the results. 



IX. RESULTS ON THE SM HIGGS BOSON 

We determine limits on the SM Higgs boson production 
cross section using a modified frequentist approach [ill 
using the distributions of the final discriminants shown 
in Fig. 131 A log-likelihood ratio (LLR) test statistic 
is formed using the Poisson probabilities for estimated 
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FIG. 3: (color online). Distribution of the BDT outputs for the (a) efifj,A, (b) efifiB, (c) efific, and (d) ee/x channels and 
distributions of the outputs of the (e,g) first and (f,h) second BDT in the HThTh and e^fj,^ channels, respectively. The data are 
compared to the sum of the expected background and to simulations of a Higgs boson signal for Mr = 125 GeV, multiplied 
by factors of 75 for the efj,fj,B and e/x/ic channels and 10 for the other channels. 



background yields, the signal acceptance, and the ob- 
served number of events for different Higgs boson mass 
hypotheses. The confidence levels are derived by in- 
tegrating the LLR distribution in pseudo-experiments 
using both the signal-plus-background (CLsb) and the 
background-only hypotheses (CLb). The excluded pro- 
duction cross section is taken to be the cross section for 
which the confidence level for signal, CL^ =CLsh/ CLt,, 
equals 0.05. The limits on the cross section for the dif- 
ferent final states are given in Table IIVI The individual 
channels have similar sensitivity, and the combined up- 
per limits only vary within about a factor of two over 
the entire mass range of 100 < Mh < 200 GeV. At 
Mh — 125 GeV the expected and observed upper Hmits 
on the cross section, expressed as a ratio relative to the 
predicted SM cross section, are 6.3 and 8.4, respectively. 
Figure S] shows the limits on the cross section and the 
LLR distributions for each channel and for the combined 
result. 



X. RESULTS ON FERMIOPHOBIC HIGGS 
BOSONS 



In addition, we set limits in a fermiophobic Higgs bo- 
son model, where the Higgs boson is assumed to couple to 
W and Z bosons with SM strengths and the Higgs boson 
couplings to fermions are zero. The gluon-gluon fusion 
Higgs production cross section is therefore small and is 
neglected. The outputs of the BDTs trained using a SM 
Higgs boson signal are used to set the limits with signal 
distributions where the gluon-gluon fusion processes and 
the H ^ T~^T~ decays have been removed and only the 
VH production mechanism is considered. The limits on 
the cross section, given as a ratio relative to the cross 
section in the fermiophobic model, are listed in Table [V] 
The ee/i, e/i/i, and e*/i^ channels have similar sensitivity 
to a fermiophobic Higgs boson, whereas the jiThTh chan- 
nel is less sensitive because H t~^t~ decays do not 
contribute. The combined expected upper limits vary 
between 2.3 and 8.5 and the observed limits between 2.4 
and 13.0, expressed as a ratio relative to the cross section 
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FIG. 4: (color online). Upper limit on the SM Higgs boson production cross section expressed as a ratio to the SM prediction 
(left column) and observed LLR (right column) as a function of Mh for the (a,b) ee/i, (c,d) e/i/i, (e,f) ^ThTh, (g,h) e*/i* 
channels, and (i,j) for all channels combined. The LLRs are shown for the background-only and the signal-plus-background 
hypotheses. The bands correspond to regions of ±1 and ±2 standard deviations (s.d.) around the median expected limit and 
around the expected median LLR for the background-only hypothesis, respectively. 
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TABLE III: Numbers of events in data, predicted background, 
and expected signal for Mh = 125 GeV for the e*^^ channel. 
The numbers are given after the initial event selection and 
after the final selection, which also requires the first BDT 
output to be > 0.3 and min{MT(e), Mt(m)} > ^ GeV. The 
numbers of events are given with their total (statistical and 
systematic) uncertainties. 



Signal 

WH 

ZH 

gg H ^ ZZ 



Signal Sum 



Background 

Z e+e" 

Z — >■ fl^ j-L^ 

Z T+r- 
Zj 

Diboson 

tt 

W+jets 
Muhijet 



Background 
Sum 



Data 



Initial Selection 



1.93 
0.32 
< 0.01 



2.25 



15.9 ±2.4 
58.5 ± 15.2 
22.0 ±6.8 

< 0.1 
36.2 ±3.6 
4.1 ± 2.1 
238.3 ± 19.0 
434.5 ± 87.0 



809 ± 93 



822 



Final Selection 



1.51 
0.23 
< 0.01 



1.74 



2.7 ±0.4 
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1.8 ±0.6 
< 0.1 

31.6 ±3.2 
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FIG. 5: (color online), (a) Upper limit on the fermiophobic 
Higgs boson production cross section expressed as a ratio to 
the prediction and (b) observed LLR as a function of Mh for 
the combined channels. Also shown are the expected LLR 
distributions for the background-only hypothesis and for the 
signal+background hypothesis, with the bands indicating ±1 
and ±2 s.d. fluctuations around the expected median LLR for 
the background-only hypothesis. 



TABLE IV: Expected and observed 95% C.L. upper limits 
on the SM Higgs boson production cross section relative to 
the SM expected value, for the ee/^, e/i/i, fj,ThTh, and 
channels separately and combined. 



eefi 



Exp Obs 



16.6 
17.4 
18.6 
17.7 
16.5 
14.1 
12.3 
11.0 
10.1 
9.4 
8.9 
9.3 
9.6 
9.6 
11.0 
11.9 
12.9 
13.6 
14.4 
15.5 
16.4 



36.1 
36.1 
34.8 
34.1 
28.6 
19.9 
17.4 
16.0 
12.6 
11.2 
11.7 
11.5 
10.7 
9.3 
10.7 
10.6 
11.3 
13.0 
13.1 
13.8 
11.6 



Exp Obs 



24.8 32.9 
23.5 24.0 
24.0 38.2 
22.3 27.1 
21.7 22.5 
17.0 22.3 
14.3 15.4 



13.1 
11.4 
11.2 
10.6 
10.8 
10.9 
10.3 
11.0 



12.4 
11.3 
11.4 
9.8 
9.0 
9.3 
8.5 
12.3 



12.7 22.4 
13.5 16.7 
14.1 19.8 
15.1 29.1 
17.1 25.6 

17.8 23.7 



fJ-ThTh 



Exp Obs 



8.2 
9.3 
10.2 
11.3 
12.7 
13.0 
13.5 
14.6 
14.1 
14.2 
16.2 
15.4 
15.4 
16.1 
16.0 
17.4 
21.1 
20.1 



10.8 
11.4 
12.3 
13.6 
17.2 
19.4 
13.3 
17.6 
20.6 
22.3 
20.1 
17.6 
22.8 
23.9 
16.2 
34.3 
40.7 
26.2 



24.0 36.7 
25.9 37.8 
27.5 33.3 



Exp Obs 



18.6 
19.3 
18.9 
17.8 
14.4 
11.8 
10.4 
8.4 
8.5 
7.7 
7.0 
7.3 
6.9 
6.6 
7.5 
7.9 
8.5 
9.9 
10.7 
12.8 
11.9 



10.4 
12.3 
13.0 
12.9 
9.8 
8.8 
7.4 
6.2 
7.2 
6.4 
6.9 
6.2 
5.9 
6.3 
7.2 
8.3 
10.4 
11.3 
14.5 
17.3 
17.3 



Combined 



Exp Obs 



6.3 
6.7 
7.4 
7.1 
7.3 
6.3 
5.9 
5.1 
4.9 
4.6 
4.3 
4.4 
4.2 
4.1 
4.5 
5.0 
5.6 
6.0 
6.7 
7.7 
7.7 



7.5 
7.2 
7.2 
10.9 
9.6 
8.4 
5.5 
4.9 
5.2 
5.1 
5.2 
4.5 
4.4 
4.6 
4.7 
7.7 
7.3 
10.6 
11.2 
12.7 
10.1 



in the fermiophobic model. Figure [5] shows the limits on 
the cross section for the production of a fermiophobic 
Higgs boson and the LLR distribution. 



XI. CONCLUSION 

We have presented the first search for the SM Higgs 
boson in multilepton eefi, e^ifi, and fir^Th final states and 
in like-charge e^/^* final states using the DO detector. 
The search is mainly sensitive to associated production 
of a, W 01 Z boson with a Higgs boson, where the Higgs 
boson decays into W'^W~ and ZZ pairs, thereby probing 
the HVV coupling directly in production and decay. We 
set limits on the cross section for a SM Higgs boson in the 
range 100 < Mh < 200 GeV with expected and observed 
upper limits of 6.3 and 8.4 at Mh = 125 GeV, expressed 
as the ratio relative to the predicted SM cross section. 
We also interpret the data in a fermiophobic Higgs boson 
model. 
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Mh 
(GeV) 




efifj, 




e± 


M± 


Combined 


Exp Obs 


Exp Obs 


Exp Obs 


Exp Obs 


Exp 


Obs 


100 


5.1 


10.4 


5.7 


7.6 


26 


27 


3.7 


2.1 


2.6 


2.5 
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4.4 
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5.1 
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17 
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18 


17 


4.1 


3.0 
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6.1 


17 


20 


4.1 


3.1 


2.8 


2.6 


140 


6.7 


8.3 


7.2 


6.8 


17 


18 


5.1 


4.1 


3.2 


2.8 


145 


7.2 


8.3 


8.1 


8.3 


16 


20 


5.3 


4.5 


3.4 


3.3 


150 


7.7 


9.8 


8.9 


7.6 


20 


28 


5.6 


5.4 


3.5 


4.0 


155 


8.4 


10.7 


9.7 


7.9 


21 


23 


6.5 


5.4 


4.0 


3.8 


160 


9.4 


10.5 


10.4 


9.0 


19 


30 


6.6 


5.7 


4.1 


4.2 


165 


9.4 


9.3 


10.0 


8.3 


22 


32 


6.6 


6.1 


4.2 


4.0 


170 


11.2 


10.8 


11.2 


12.1 


23 


30 


7.3 


7.0 


4.7 


4.9 


175 


12.1 


10.6 


12.7 


21.9 


25 


29 


7.9 


8.2 


5.1 


6.5 


180 


13.5 


11.4 


14.4 


16.4 


29 


33 


8.5 


10.1 


5.8 


6.5 


185 


14.8 


13.7 


15.3 


20.6 


31 


45 


10.0 


11.4 


6.4 


8.5 


190 


16.3 


13.5 


17.0 


30.5 


35 


50 


11.0 


14.4 


7.2 


10.8 


195 


17.8 


15.5 


19.7 


29.1 


40 


62 


12.9 


17.3 


8.2 


13.0 
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19.0 


12.0 


20.6 


23.9 


43 


47 


12.1 


17.3 


8.5 


9.9 
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TABLE V: Expected and observed 95% C.L. upper limits on 
the fermiophobic Higgs boson production cross section rela- 
tive to the expected cross section in the fermiophobic model, 
for the eefi, e/i/i, fiThTh, and e^^t* channels separately and 
combined. 
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